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Abstract. In this work we develop a lepto-hadronic model for the elec- 
tromagnetic radiation from jets in microquasars with low-mass compan- 
ion stars. We present general results as well as applications to some spe- 
cific systems, and carefully analyze the predictions of the model in the 
gamma-ray band. The results will be directly tested in the near future 
with the present and forthcoming space-borne and terrestrial gamma-ray 
telescopes. 

Resumen. En este trabajo se desarrolla un modelo lepto-hadronico pa- 
ra la radiacion electromagnetica de jets en microcuasares con estrellas 
compaheras de baja masa. Se presentan resultados generales y aplicacio- 
nes a sistemas especfiicos, y se analizan en detalle las predicciones del mo- 
delo en la banda de rayos gamma. Los resultados podran ser directamente 
contrastados en el futuro cercano con las observaciones de telescopios de 
rayos gamma espaciales y terrestres de presente y futura generacion. 

1. Introduction 

An outstanding feature of astrophysical accreting sources at all scales is the pro- 
duction of jets - collimated, bipolar, extended flows of matter and electromag- 
netic field ejected from the surroundings of a rotating object. Jets are launched 
from accreting supermassive and stellar-mass black holes, neutron stars, white 
dwarfs, and protostars (Romero et al. 2010). X-ray binaries (XRBs) with jets 
are called microquasars (Mirabel et al. 1992). Microquasars are formed by a 
non-collapsed star and a stellar-mass compact object, that may be a neutron 
star or a black hole. The compact object accretes matter lost by the companion 
star. A fraction of this matter is ejected from the system as two collimated jets. 

According to the mass of the donor star, microquasars (and all XRBs) are 
classified into low-mass and high-mass. In high-mass microquasars (HMMQs) 
the donor star is an O, B, or Wolf Rayet star of mass ~ 8 — 2OM . These stars 
lose mass mainly through strong winds. Donor stars in low-mass microquasars 
(LMMQs) have masses < 2M . They are old stars of spectral type B or later, 
that transfer mass to the compact object through the overflow of their Roche 
lobe. The position of the known galactic microquasars (~ 20 sources) is shown 
in Figure 1. High- mass microquasars trace the star forming regions in the spi- 
ral arms of the galaxy (Bodaghee et al. 2012, Coleiro & Chaty 2011). This is 
expected, since the companion star is relatively young (< 10 7 yr) and should 
not have departed significantly from its birthplace. Low-mass microquasars are 
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concentrated towards the center of the galaxy, especially in the bulge. These 
systems are old (~ 10 9 yr), and some have migrated from the galactic plane 
towards higher latitudes (e.g. Mirabel et al. 2001). 
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Figure 1. Spatial distribution in galactic coordinates of microquasars in 
the Milky Way. Data of 2010 collected by Dr. S. Chaty, available online at 
http : //www. aim. univ-paris7 . f r/CHATY/Microquasars/microquasars .html. 
The nature of some of the systems (e.g. LS 5039 and LS I +61° 303) is still 
disputed, but they have been included in the figure for historical reasons. 

Black hole X-ray binaries go through different spectral states, classified 
according to the timing and spectral characteristics of the X-ray emission (e.g. 
McClintock & Remillard 2006, Belloni et al. 2011). In microquasars, steady jets 
are observed in the low-hard state (and probably also in quiescence), whereas 
discrete jets are ejected during the transition between states. Steady jets have 
typical luminosities ~ i0 36-37 er g g -i^ anc j are ^idly relativistic with bulk 
Lorentz factors ~ 1.5 — 2. Discrete ejections may be much faster: bulk velocities 
close to the speed of light are inferred from the apparent superluminal motion 
of some blobs (e.g. Mirabel & Rodriguez 1994). 

The electromagnetic emission of steady jets extends from radio wavelengths 
up to a turnover at infrared/optical frequencies. The spectrum is clearly non- 
thermal and very well explained as synchrotron radiation of a population of 
relativistic electrons. Although the - also non-thermal - hard X-ray emission in 
XRBs is generally attributed to the corona, in some microquasars a correlation 
between radio and X-rays is observed during the low-hard state (Corbel et al. 
2003, Gallo et al. 2003). This suggests that the radiation in both bands orig- 
inates in the same region (the jets) and by the same process (electron syn- 
chrotron). The HMMQs Cygnus X-1 (McConnell et al. 2000, Albert et al. 2007, 
Sabatini et al. 2010) and Cygnus X-3 (Tavani et al. 2009, Abdo et al. 2009) have 
been detected at gamma-ray energies. This emission is thought to be produced 
also in the jets through the interaction of relativistic particles with the wind 
and the radiation field of the companion star (e.g. Bosch-Ramon et al. 2008, 
Araudo et al. 2010, Romero et al. 2010). No LMMQ (or low-mass XRB) has 
been detected in gamma rays yet. In these systems the old and dim companion 
star cannot provide enough targets to the jets. Then, if LMMQs were to be 
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high energy sources, we expect the gamma rays to be created in the outflows by 
interaction of relativistic particles with the magnetic field, matter, and radiation 
inside the jets. 1 

In this work we develop a model for the electromagnetic radiation of jets 
in microquasars that generalizes and improves those existing in the literature. 
We seek to obtain a better understanding of the physical conditions in the jets 
through the comparison of our results with observational data. We focus in 
particular on models for jets in low-mass microquasars, seeking to assess their 
detectability at high (~ GeV) and very high (~ TeV) energies with the instru- 
ments available now or in the near future. In the following sections we outline the 
jet model and present some general results and applications to specific sources. 

2. Jet model 

The jet model is described in detail in Romero & Vila (2008), Vila & Romero 
(2010), Vila et al. (2012), and Vila (2012). 2 Here we only review it briefly. 
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Figure 2. Left: sketch of a microquasar. Right: a detail of the jet. Some 
relevant geometrical parameters are indicated. 

A sketch of the system under study is shown in Figure 2. We assume that 
the jets are launched at a distance zo = 50i? grav 3 from the compact object 
with an initial radius tq = O.lzcb and expand conically while advancing with 
a constant bulk Lorentz factor Rjet- The axis of the jet makes an angle #j e t 
with the line of sight. The power of each jet is taken to be a fraction of the 
accretion power, L- ]et = q ]et L &ccr . The hypothesis of equipartition between the 
magnetic energy density and the kinetic energy density of the outflow allows 



The corona might also emit gamma rays, see for example Vieyro & Romero (2012). 

2 The complete manuscript of the PhD Thesis (Vila 2012) is available for download at 
http: //f caglp .unlp . edu. ar/~gvila/. 

3 The gravitational radius of a black hole of mass M is -R gra v = GM/c 2 . 



4 



G. S. Vila 



to estimate the magnetic field at the base of the jet. We obtain typical values 
So = B(zq) ~ 10 5 ~ 6 G. The field strength decreases as the jet expands; we adopt 
the prescription B(z) = BQ(zo/z) m with m > 1 to parameterize its evolution. 

The composition of the jets is unknown, although the existence of acceler- 
ated electrons (or electron-positron pairs) is firmly inferred from the observed 
synchrotron radio spectrum. In at least one microquasar (SS433, Migliari et al. 
2002) there is evidence of the presence of nuclei as well. We assume that the 
outflows contain electrons and protons, and that both species are accelerated 
up to relativistic energies with the same efficiency. In our model a fraction 
-^rei = QjetLjet of the jet power is transferred to relativistic particles; we fix 
Qrei = 0.1. The power is shared between protons and electrons, so L rc \ = L p + L e 
with L p = aL e . We consider values of a > 1 (proton-dominated jets) as well as 
the case a = 1 (equipartition of energy between both species). 

Particles are accelerated in the region of the jet z acc < z < z max . We do 
not investigate the acceleration mechanism, but assume that the injection spec- 
trum follows a power-law in energy, Q CjP oc E~^ (in units of s -1 cm -3 erg" 1 ) 
with 1.5 < a < 2.4. This is consistent with acceleration by diffusion across 
shock fronts (also known as first order Fermi process, see e.g. Drury 1983). 
Particles gain energy but simultaneously cool. The maximum energy they can 
reach is calculated demanding that the acceleration rate equals the total energy 
loss rate, t~,} c (E mauX ) = t~^ ol (E max ). In the cooling rate we include adiabatic 
energy losses and several radiative processes: synchrotron radiation, relativis- 
tic Bremsstrahlung, inverse Compton interaction, proton-proton, and proton- 
photon inelastic collisions. Protons achieve maximum energies much larger than 
electrons. 

Proton-proton and proton-photon inelastic collisions produce pions. Neu- 
tral pions decay into two gamma rays, whereas charged pions decay into muons 
and muon neutrinos. The subsequent decay of muons injects electrons (or 
positrons) and more neutrinos. A second channel of proton-photon interaction 
directly creates electron-positron pairs, too. We call the charged products of 
hadronic interactions "secondary" particles to distinguish them from "primary" 
protons and electrons - those that are directly accelerated in the jet as discussed 
in the previous paragraph. 1 The mean lifetime of charged pions and muons is 
very short, so the possibility that they cool significantly before decaying is usu- 
ally disregarded. If these particles are created energetic enough, however, the 
synchrotron cooling time in the strong magnetic field of the jet may be shorter 
than their mean lifetime. We consider the cooling of all (stable and unstable) 
secondary particles when computing their steady-state distribution (see below) 
and calculate their contribution to the radiative spectrum of the jet through the 
same radiative mechanism than for primary particles. 

For all (primary and secondary) particle species the steady-state energy 
distribution N(E,z) (in units of cm -3 erg -1 ) is calculated solving a transport 
equation that takes into account injection, cooling, convection or escape, and 
particle decay. We work with two versions of the transport equation that repre- 
sent different assumptions about the conditions in the acceleration region. The 



4 We also include among the secondary particles the pairs injected through two-photon annihi- 
lation, see Section 3. 
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initial models are "one- zone". In these the acceleration region is compact and 
homogeneous, and placed either at the base of the jet or at some other posi- 
tion z acc where the kinetic energy density dominates over the magnetic energy 
density. Magnetohydrodynamical models predict that this condition favors the 
formation of shock waves. In one-zone models the removal of particles from the 
acceleration region is accounted for through an escape term in the transport 
equation. Later, we study the injection and cooling of relativistic particles in a 
spatially extended, inhomogeneous region of the jet. In these models a connec- 
tive term with a spatial derivative is added to the transport equation and the 
escape term is removed. 

3. General results 

Figure 3 shows the spectral energy distributions (SEDs) obtained for two rep- 
resentative sets of values of the parameters of the one-zone jet model; see 
Romero & Vila (2008) for details. The left panel corresponds to a model with 
a = 1 and the right panel to a proton-dominated jet with a = 10 3 . In both 
cases the emission from radio to hard X-rays is synchrotron radiation of rela- 
tivistic electrons and the gamma rays are of hadronic origin (decay of tt° and 
synchrotron radiation of created via jry collisions). The most relevant pho- 
ton target field for jrf interactions is the synchrotron field of electrons. Thus, 
in the model with a = 1 we obtain appreciable very high-energy gamma-ray 
luminosity at the level of ~ 10 35 erg s" 1 . In the proton-dominated model, on 
the other hand, the low luminosity of the electron synchrotron component con- 
sequently quenches the very high-energy emission - this even though 10 3 times 
more power is injected in relativistic protons than in electrons. Notice, however, 
the prominent proton synchrotron peak at ~ 1 GeV. 




Log 10 [£, /eV] Log l0 [E., /eV] 

Figure 3. SEDs calculated in a one-zone jet model with a = 1 and a = 2.2 
(left), and a = 10 3 and a = 1.5 (right). 

The high-energy radiation may be absorbed before escaping the jet. We 
consider only one absorption process, the annihilation of two photons into an 
electron-positron pair. Figure 4 shows the SEDs corrected by absorption. Gamma 
rays annihilate efficiently against IR and UV photons, that once again are pro- 
vided by the electron synchrotron radiation field. All emission above ~ 1 GeV 
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is therefore suppressed in the model with a = 1. A source like this would not be 
detectable by Cherenkov telescopes at very high energies, although it might be 
marginally detected at high energies with a soft spectrum. The SED is unmodi- 
fied in the case a = 10 3 . Our model predicts that a source of these characteristics 
might be detectable at GeV energies and possibly at TeV energies depending on 
the distance. 

From these general examples it is clear that considering purely leptonic or 
purely hadronic models is hardly meaningful, since the efficiency of hadronic 
processes to produce gamma rays and the opacity to gamma-ray propagation 
strongly depend on the radiation field of electrons at lower energies. A second 
point is worth remarking: a neutrino luminosity approximately of the same order 
than the gamma-ray luminosity is injected through jry interactions. Neutrinos, 
however, freely escape the jet. According to our model, then, sources strongly 
absorbed in gamma rays should be the most luminous neutrino emitters - no 
correlation between gamma rays and neutrinos should be a priori expected. 




Log 10 [E, /eV] Log,„ [E y /eV] 



Figure 4. The SEDs in Figure 3 modified by absorption. The sensitivities of 
Fermi (5a, one- year sky survey exposure), HESS (5a, 50 h exposure), MAGIC 
II (50 h exposure), and the predicted for CTA (50 h exposure) are indicated. 
In the right panel, the sensitivity curves of HESS and MAGIC II are plotted 
for a source at 2 kpc (higher sensitivity) and 6 kpc (lower sensitivity). 



4. Some specific applications 

We applied the model to the study of two black hole LMMQs, the sources 
GX 339-4 (Vila & Romero 2010) and XTE J1118+480 (Vila et al. 2012). Both 
are X-ray transients: they spend years in quiescence until due to a sudden 
increase in the mass accretion rate they enter in outburst for some months. Five 
outbursts have been observed in GX 339-4 and two in XTE J1118+480. In all 
occasions they were extensively and simultaneously (or quasi-simultaneously) 
monitored in several bands from radio to X-rays. The X-ray spectrum was 
consistent with the sources being in the low-hard state. A radio jet was imaged 
in GX 339-4 (Gallo et al. 2004) but not yet in XTE Jll 18+480, although its 
presence is inferred from the characteristics of the radio emission. 
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Figure 5 shows one-zone model fits to broadband observations of GX 339-4 
during the outburst of 1997 (left panel) and a low luminosity phase of the out- 
burst of 1999 (right panel). For these applications we located the base of the 
acceleration region far from the base of the jet, typically z acc ~ 10 2 — 10 4 i? grav . 
The emission from radio to X-rays is synchrotron radiation of primary electrons 
and, in the fit to the data from 1997, secondary electron-positron pairs created 
by photon-photon annihilation. In this case there is also significant synchrotron 
emission of pions and muons, although not dominant in any energy range. In 
the fit to the low- luminosity SED from 1999 a data point in the optical was 
included in the fit. The rising shape of the spectrum at these energies cannot be 
explained as emission from the jet alone, and possibly originates in an accretion 
disk. Both theoretical SEDs clearly differ in the predictions for the gamma- 
ray band. Whereas in a high X-ray luminosity outburst like that of 1997 the 
source might be detectable at high and very high energies, the model predicts 
no detection during a low-luminosity state like the observed in 1999. A word of 
caution is, however, due here: the same data may be fit with equal success us- 
ing different sets of values of the model parameters, see Vila & Romero (2010). 
This degeneracy is inherent to the modeling. We expect it to be removed in the 
near future when gamma rays observations of the source become available, since 
different fits predict different spectral shapes at high and very high energies. 




Figure 5. Model fits to broadband observations of GX 339-4 during the 
outbursts of 1997 (left panel) and 1999 (right panel). Data collected from 
various authors, see Vila & Romero (2010) for details. 



We performed fits to broadband observations of XTE J1118+480 during 
the outbursts of 2000 and 2005; the results are plotted in Figure 6. The SEDs 
were calculated applying the inhomogeneous jet model. A thermal component 
peaking at ~ 20 — 40 eV is detected in this microquasar, so we added a sim- 
ple geometrically thin, optically thick accretion disk model (Shakura & Sunyaev 
1973) to our calculations. According to our modeling, there are no great dif- 
ferences in the physical conditions in the jets during the two outbursts. The 
radio and X-ray emission is fitted by the synchrotron spectrum of primary elec- 
trons, plus some contribution at low energies of secondary pairs created through 
photon-photon annihilation; the emission of other species of charged secondaries 
is negligible. The IR-optical-UV range has significant contribution from the ac- 
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cretion disk. The IC scattering off the jet photon field by primary electrons 
contributes in a narrow energy range about ~ 10 GeV in the case of the 2000 
outburst. The SED above ~ 1 GeV is completely dominated by gamma rays 
from the decay of neutral pions created in proton-proton collisions. The photon 
field of the disk is the main source of opacity to photon escape. The optical 
depth is large only near the base of the acceleration region, z ~ -^acc ^ io 8 cm; 
gamma rays with energies 10 GeV < -E7 7 < 1 TeV are mostly absorbed in this 
zone. The total luminosities are, nevertheless, unmodified by absorption. The 
reason is that there are many high-energy protons that produce gamma rays 
through proton-proton collisions outside the acceleration region. This radiation 
escapes unabsorbed since the density of disk photons is low at high z. According 
to our results, a future outburst of the source with emission levels comparable 
to those of 2000 and 2005, would be detectable in gamma rays by ground-based 
observatories like MAGIC II and CTA. In the context of the model presented 
here, observations at very high energies would help to constrain the hadronic 
content of the jets, since above ~ 100 GeV the predicted emission is completely 
due to proton- proton interactions. 
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Figure 6. Model fits to broadband observations of XTE J1118+480 dur- 
ing the outbursts of 2000 (left panel) and 2005 (right panel). Data from 
McClintock et al. (2001) and Maitra ct al. (2009). 



5. Final comments and perspectives 

Along this work we developed a lepto-hadronic model for the radiation from jets 
in microquasars. In its current version the model is suited to the study of low- 
mass microquasars, but there are a number of specific improvements that would 
allow a broader application. We have attempted to add as much detail and 
self-consistency to the model as possible. It depends on a number of parameters 
for which we took, when available, estimates inferred from observations. The 
values of other parameters (such as the position of the acceleration region and 
the spectral index of the particle injection function) were chosen to account, in 
an affective manner, for the main constraints imposed by the physics of particle 
acceleration and the dynamics of the outflow. Our modeling is, nevertheless, by 
no means free of limitations. The most important one is, perhaps, related to 
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the magnetic field. We have made simple and reasonable assumptions on this 
point, but they might require further refinement in the light of future insights 
and observational data. 

We have shown that the spectral energy distributions from microquasar 
jets might be complex and take a variety of shapes depending on the condi- 
tions in the source; we have tried to cover a large number of scenarios within a 
physically meaningful range for the values of the model parameters. Relativis- 
tic jets from low-mass X-ray binaries with a content of non-thermal electrons 
and protons radiate along the whole electromagnetic spectrum. From radio to 
X-rays the emission is of leptonic origin, predominantly due to synchrotron ra- 
diation. In proton-dominated models (a > 1), the spectrum above ~ 1 GeV 
is of hadronic origin (synchrotron, proton-proton, and proton-photon inelastic 
collisions). Under some particular conditions the high and very high energy 
gamma-ray emission might reach levels detectable with presently operative in- 
struments. Proton-proton and proton-photon inelastic collisions inject charged 
pions and muons. For values of the magnetic field as those we adopted the cool- 
ing of these particles before decay is not negligible. Although radiatively not 
relevant, the cooling of pions and muons seriously affect the neutrino emissiv- 
ity from the jets, see Reynoso & Romero (2009). Gamma rays can escape the 
source without significant absorption if the emission region is located in a zone 
of the jet where the radiation field has low density. In one-zone models this 
is possible if the relativistic particles are injected at large distances from the 
black hole. In models for extended jets absorption can be avoided even if the 
acceleration region is relatively near the jet base: plenty of energetic protons 
leave the acceleration region and inject gamma rays in zones with low internal 
and external photon density. 

We performed fits to observations of very well studied low-mass micro- 
quasars, GX 339-4 and XTE J1118+480. We make predictions for the high 
and very high energy gamma-ray spectrum during outbursts, a question not 
addressed in previous works about these sources. The model can satisfactorily 
reproduce the data, although sometimes with more than one set of parameters. 
Some of our models predict detectable gamma-ray emission for these sources in 
outburst. 

To date no low-mass X-ray binaries have been detected at high or very 
high energies, although there are some candidates. Their observation is further 
complicated because in general they are transient sources. We expect that this 
situation changes in the very near future. The detection (or not) of low-mass 
XRBs at high and very high energies will provide very valuable information. 
The most favorable situation would be, undoubtedly, to have at our disposal 
simultaneous observations in X-rays and high and very high energy gamma rays. 
Such simultaneous spectral coverage is nowadays possible. Together with the 
much improved quality of the data, it will allow to remove part of the inherent 
degeneracy of the modeling. In this context, and in spite of its limitations, 
the type of models developed in this thesis are timely. We expect that, when 
confronted with observations, they result adequate to reproduce the radiative 
spectrum from microquasar jets and contribute to a better understanding of 
these objects. 
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